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Abstract. Heteroclinic cycles involving two saddle-foci, where the saddle-foci share both in¬ 
variant manifolds, occur persistently in some symmetric differential equations on the 3-dimensional 
sphere. We analyse the dynamics around this type of cycle in the case when trajectories near 
the two equilibria turn in the same direction around a 1-dimensional connection — the saddle- 
foci have the same chirality. When part of the symmetry is broken, the 2-dimensional invariant 
manifolds intersect transversely creating a heteroclinic network of Bykov cycles. 

We show that the proximity of symmetry creates heteroclinic tangencies that coexist with 
hyperbolic dynamics. There are n-pulse heteroclinic tangencies — trajectories that follow the 
original cycle n times around before they arrive at the other node. Each n-pulse heteroclinic 
tangency is accumulated by a sequence of (n-l-l)-pulse ones. This coexists with the suspension of 
horseshoes defined on an infinite set of disjoint strips, where the first return map is hyperbolic. 
We also show how, as the system approaches full symmetry, the suspended horseshoes are 
destroyed, creating regions with infinitely many attracting periodic solutions. 


1. Introduction 

A Bykov cycle is a heteroclinic cycle between two hyperbolic saddle-foci of different Morse 
index, where one of the connections is transverse and the other is structurally unstable — see 
Figure [T] There are two types of Bykov cycle, depending on the way the flow turns around 
the two saddle-foci, that determine the chirality of the cycle. Here we study the non-wandering 
dynamics in the neighbourhood of a Bykov cycle where the two nodes have the same chirality. 
This is also studied in [M], and the case of different chirality is discussed in [36]. A simplified 
version of the arguments presented here appears in [35]. 

1.1. The object of study. Our starting point is a fully (Z 2 x Z 2 )-symmetric differential equa¬ 
tion X = fo{x) in the three-dimensional sphere with two saddle-foci that share all the invariant 
manifolds, of dimensions one and two, both contained in flow-invariant submanifolds that come 
from the symmetry. This forms an attracting heteroclinic network with a non-empty basin of 
attraction . We study the global transition of the dynamics from this fully symmetric system 
X = fo{x) to a perturbed system x = fx{x), for a smooth one-parameter family that breaks part 
of the symmetry of the system. For small perturbations the set is still positively invariant. 

When A 7 ^ 0, the one-dimensional connection persists, due to the remaining symmetry, and the 
two dimensional invariant manifolds intersect transversely, because of the symmetry breaking. 
This gives rise to a network that consists of a union of Bykov cycles, contained in V^. 
For partial symmetry-breaking perturbations of /o, we are interested in the dynamics in the 
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Figure 1. A Bykov cycle with nodes of the same chirality. There are two possibilities 
for the geometry of the flow around a Bykov cycle depending on the direction trajectories 
turn around the connection [v — >■ w]. We assume here that the nodes have the same 
chirality: trajectories turn in the same direction around the connection. When the 
endpoints of a nearby trajectory are joined, the closed curve is always linked to the cycle. 


maximal invariant set contained in V^. It contains, but does not coincide with, the suspension 
of horseshoes accumulating on described in [3l EU [Ml SS]. Here, we show that close to the 
fully symmetric case it contains infinitely many heteroclinic tangencies. Under an additional 
assumption, we show that Vq contains attracting limit cycles with long periods, coexisting with 
sets with positive entropy. 

1.2. History. Homoclinic and heteroclinic bifurcations constitute the core of our understanding 
of complicated recurrent behaviour in dynamical systems. This starts with Poincare on the late 
19*^ century, with major subsequent contributions by the schools of Andronov, Shilnikov, Smale 
and Palis. These results rely on a combination of analytical and geometrical tools used to 
understand the qualitative behaviour of the dynamics. 

Heteroclinic cycles and networks are flow-invariant sets that can occur robustly in dynamical 
systems with symmetry, and are frequently associated with intermittent behaviour. The rigorous 
analysis of the dynamics associated to the structure of the nonwandering sets close to heteroclinic 
networks is still a challenge. We refer to [28] for an overview of heteroclinic bifurcations and for 
details on the dynamics near different kinds of heteroclinic cycles and networks. 

Bykov cycles have been found analytically in the Lorenz model in [U |M] and the nearby 
dynamics was studied by Bykov in in miiniiii]. The point in parameter space where this cycle 
occurs is called a T-point in [T8|. Recently, there has been a renewal of interest in this type 
of heteroclinic bifurcation in the reversible usmsiiMi, equivariant laiMiiiHi and conservative 
cases [Tj. See also m- 

The transverse intersection of the two-dimensional invariant manifolds of the two equilibria 
implies that the set of trajectories that remain for all time in a small neighbourhood of the 
Bykov cycle contains a locally-maximal hyperbolic set admitting a complete description in terms 
of symbolic dynamics, reminiscent of the results of Shilnikov [52] . An obstacle to the global 
symbolic description of these trajectories is the existence of tangencies that lead to the birth of 
stable periodic solutions, as described for the homoclinic case in ladziiiniiiiiiM]. 

All dynamical models with quasi-stochastic attractors were found, either analytically or by 
computer simulations, to have tangencies of invariant manifolds laiiiiiig. As a rule, the sinks 
in a quasi-stochastic attractor have very long periods and narrow basins of attraction, and they 
are hard to observe in applied problems because of the presence of noise |24] . 

Motivated by the analysis of Lamb et al m in the context of the Michelson system, Bykov 
cycles have been considered by Knobloch et al EH. Using the Lin’s method, the authors extend 
the analysis to cycles in spaces of arbitrary dimension, while restricting it to trajectories that 
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remain for all time inside a small tubular neighbourhood of the cycle. We also refer the reader to 
[32], where the authors consider non-elementary T-points in reversible differential equations. The 
leading eigenvalues at the two equilibria are real and the two-dimensional invariant manifolds 
meet tangentially. They found chaos in the unfolding of this T-point and bifurcations of periodic 
solutions in the process of annihilation of the shift dynamics. 

Bykov cycles appear in many applications like the Kuramoto-Sivashinsky systems ffSlET], 
magnetoconvection [461150] and travelling waves in reaction-diffusion dynamics [61126]. 

1.3. Chirality. For a Bykov cycle in a 3-dimensional manifold, we say that the nodes have 
the same chirality if trajectories turn in the same direction around the common 1-dimensional 
invariant manifold. If they turn in opposite directions we say the nodes have different chirality. 
A more formal definition, using links, will be given in Section 12.21 below, showing this to be a 
global topological invariant of the connection that is well defined only in a 3-dimensional ambient 
space. 

These cycles have been studied by different authors who were not aware of the chirality, 
ignoring what looked like a very small and unimportant choice in local coordinates. For instance, 
Bykov in mm addresses the case of different chirality without mentioning it explicitly — see 
a discussion in Section 7 of [36] . Cycles with the same chirality are treated in [SUM] and they 
occur naturally in reversible differential equations [32] . Dynamical features that are irrespective 
of chirality are described in m- 

Arbitrarily close to Bykov cycles of any chirality there are suspended horseshoes and multi¬ 
pulse heteroclinic cycles — see Theorem [2| below. Around Bykov cycles where the nodes have 
different chirality, heteroclinic tangencies occur generically in trajectories that remain close to the 
cycle for all time, as shown in [36] . This is not the case when the nodes have the same chirality, 
but we show here that heteroclinic tangencies appear when the equations are approaching a 
more symmetric one. 

1.4. Symmetry. Heteroclinic cycles involving equilibria are not a generic feature in differential 
equations, but they can be structurally stable in systems which are equivariant under the action 
of a symmetry group, due to the existence of flow-invariant subspaces [26]. Thus, perturbations 
that preserve the symmetry will not destroy the cycle. Explicit examples of equivariant vector 
fields for which such cycles may be found are reported in [lElEniEellMlIiaile]. Symmetry, 
exact or approximate, plays an essential role in the analysis of nonlinear physical phenomena 
[6] da He]. It is often incorporated in models either because it occurs approximately in the 
phenomena being modelled, or because its presence simplifies the analysis. Since reality has not 
perfect symmetry, it is desirable to understand the dynamics that are being created under small 
symmetry breaking perturbations. 

Symmetry plays two roles here. First, it creates flow-invariant subspaces where non-transverse 
heteroclinic connections are persistent, and hence Bykov cycles are robust in this context. Sec¬ 
ond, we use the proximity of the fully symmetric case to capture more global dynamics. Sym¬ 
metry constrains the geometry of the invariant manifolds of the saddle-foci and allows us some 
control of their relative positions, and we find infinitely many heteroclinic tangencies correspond¬ 
ing to trajectories that make an excursion away from the original cycle. For Bykov cycles with 
the same chirality, tangencies only take place near symmetry. 

In the analysis of the annihilation of hyperbolic horseshoes associated to tangencies, on the 
one hand, the symmetry adds complexity to the problem, because the analysis is not so standard 
as in [431154] . On the other hand, symmetry simplifies the analytic expression of the return map. 
It is clear that, for Bykov cycles of the same chirality, the annihilation of hyperbolic horseshoes 
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associated to tangencies only takes place near symmetry. In the fully asymmetric case, the 
general study seems to be analitically untreatable. 

Many questions remain for future work. One obvious question is to assume different chirality 
of the nodes. Due to the infinite number of reversals described in [36] , other types of bifurcations 
may occur. The second question is to observe whether the non-wandering set may be reduced 
to a homoclinic class. 

1.5. This article. We study the dynamics arising near a symmetric differential equation with 
a specific type of heteroclinic network. We show that when part of the symmetry is broken, 
the dynamics undergoes a global transition from hyperbolicity coexisting with infinitely many 
sinks, to the emergence of regular dynamics. We discuss the global bifurcations that occur as 
a parameter A is used to break part of the symmetry. We complete our results by reducing 
our problem to a symmetric version of the structure of Pahs and Takens’ result |43l §3] on 
homoclinic bifurcations. Being close to symmetry adds complexity to the dynamics. Chirality 
is an essential information in this problem. 

This article is organised as follows. In Section [21 after some basic definitions, we describe 
precisely the object of study and we review some of our recent results related to it. In Section [3| 
we state the main results of the present article. The coordinates and other notation used in 
the rest of the article are presented in Section where we also obtain a geometrical description 
of the way the flow transforms a curve of initial conditions lying across the stable manifold 
of an equilibrium. In Section [5l we prove that there is a sequence of parameter values A* 
accumulating on 0 such that the associated flow has heteroclinic tangencies. In Section [6l we 
discuss the geometric constructions that describe the global dynamics near a Bykov cycle as the 
parameter varies. We also describe the limit set that contains nontrivial hyperbolic subsets and 
we explain how the horseshoes disappear as the system regains full symmetry. We show that 
under an additional condition this creates infinitely many attracting periodic solutions. 

2. The object of study and preliminary results 

In the present section, after some preliminary definitions, we state the hypotheses for the 
system under study together with an overview of results obtained in |34j . emphasizing those 
that will be used to explain the loss of hyperbolicity of the suspended horseshoes and the 
emergence of heteroclinic tangencies near the cycle. 

2.1. Definitions. Let / be a vector held on the unit three-sphere S^, r > 3, with how given 
by the unique solution x{t) = ip{t,xo) G of 

(2.1) X = f\{x) and x(0) = xq- 

where r > 3 and A € R. Suppose that v and w are two hyperbolic saddle-foci of (12.11) with 
different Morse indices (dimension of the unstable manifold), say 1 and 2. There is a heteroclinic 
cycle associated to {v, w} if 

W“(v) n W%w) / 0 and IT“(w) n IT*(v) / 0 

where IT^(*) and IT“(a) refer to the stable and unstable manifolds of the hyperbolic saddle a, 
respectively. The terminology [v —>■ w] or [w —)• v] denotes a solution contained in IL“(v) n 
IT^(w) or IT“(v) n IT®(w), respectively. A heteroclinic network is a hnite connected union of 
heteroclinic cycles. 

For A = 0, there is a 1-dimensional trajectory in IT“(v) n VF^(w) and a 2-dimensional con¬ 
nected how-invariant manifold contained in IT“(w)nIT^(v), meaning that there are a continuum 


GLOBAL BIFURCATIONS CLOSE TO SYMMETRY March 7, 2016 


5 


of solutions connecting w and v. For A 7 ^ 0, the one-dimensional manifolds of the equilibria co¬ 
incide and the two-dimensional invariant manifolds have a transverse intersection. This second 
situation is what we call a Bykov cycle. 

These objects are known to exist in several settings and are structurally stable within certain 
classes of G-equivariant systems, where G C 0(n) is a compact Lie group. Here we consider 
differential equations ( 12 .ip with the equivariance condition: 

fxijx) = 'yfxix), for all 7 G G, A G R. 

An isotropy subgroup of G is a set G = {7 G G : 7 X = x} for some x in phase space; we write 
Fix(G) for the vector subspace of points that are fixed by the elements of G. For G-equivariant 
differential equations each subspace Fix(G) is flow-invariant. The group theoretical methods 
developed in m are a powerful tool for the analysis of systems with symmetry. 

Suppose there is a cross-section S to the flow of ( 12 . 11 ) . such that S contains a compact 
invariant set A where the first return map is well defined and conjugate to a full shift on a 
countable alphabet, we call the flow-invariant set A = t G R, g G A} a suspended 

horseshoe. 

2.2. The organising centre. The starting point of the analysis is a differential equation x = 
fo{x) on the unit sphere = {X = (xi,X 2 ,X 3 ,X 4 ) G R^ : ||A|| = 1} where /o : TS^ is a 

vector held with the following properties: 

(PI) The vector held /o is equivariant under the action of G = Z 2 © Z 2 on induced by the 
linear maps on R'^: 

7 i(xi,X2 ,X3,X4) = (-Xi,-X2,X3,X4) and 72(2:1, X2,X3,X4) = (xi,X2,-X3,X4). 

(P2) The set Fix(Z 2 ©Z 2 ) = {x G : 71 X = 72 X = x} consists of two equilibria v = (0,0,0,1) 
and w = (0,0,0, —1) that are hyperbolic saddle-foci, where: 

• the eigenvalues of d/o(v) are —C^ ± a^i and Flv with 7 ^ 0 , Cv > > 0 ; 

• the eigenvalues of dfo{w) are ± and —Gw with Ow 7 ^ 0 , Gw > E^ > 0 . 

(P3) The how-invariant circle Fix(( 7 i)) = {x G : 71 X = x} consists of the two equilibria v 

and w, a source and a sink, respectively, and two heteroclinic trajectories [v —>■ w]. 

(P4) The /o-invariant sphere Fix(( 72 )) = {x G : 72 X = x} consists of the two equilibria v 
and w, and a two-dimensional heteroclinic connection from w to v. Together with the 
connections in (Pl3|) this forms a heteroclinic network that we denote by 

Given two small open neighbourhoods V and VF of v and w respectively, consider a piece of 
trajectory ip that starts at dV, goes into V and then goes once from V to W, enters W and 
ends at dW. Joining the starting point of to its end point by a line segment, one obtains a 
closed curve, the loop of ip. For almost all starting positions in dV, the loop of ip does not meet 
the network If there are arbitrarily small neighbourhoods V and W for which the loop of 
every trajectory is linked to Sq, we say that the nodes have the same chirality as illustrated in 
Figure [TJ This means that near v and w, all trajectories turn in the same direction around the 
one-dimensional connections [v — >■ w]. This is our last hypothesis on /q: 

(P5) The saddle-foci v and w have the same chirality. 

Condition (F[5|) means that the curve ip and the cycle S® cannot be separated by an isotopy. 
This property is persistent under small smooth perturbations of the vector field that preserve 
the one-dimensional connection. An explicit example of a family of differential equations where 
this assumption is valid has been constructed in [l9]. The rigorous analysis of a case where 
property (F|5|) does not hold has been done by the authors in [36] . 
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2.3. The heteroclinic network of the organising centre. The heteroclinic connections in 

the network are contained in fixed point subspaces satisfying the hypothesis (HI) of Krupa 
and Melbourne [33]. Since the inequality C^C-w > holds, the stability criterion [33] may 

be applied to and we have: 

Lemma 1. Under conditions m-m the heteroclinic network is asymptotieally stable. 

As a consequence of Lemma [1] there exists an open neighbourhood of the network such 
that every trajectory starting in remains in it for all positive time and is forward asymptotic 
to the network. The neighbourhood may be taken to have its boundary transverse to the vector 
field /q. The flow associated to any C^-perturbation of /o that breaks the one-dimensional 
connection should have some attracting feature. 

2.4. Breaking the Z 2 (( 7 i))-symmetry. When the symmetry Z 2 (( 7 i)) is broken, the two one¬ 
dimensional heteroclinic connections are destroyed and the cycle disappears. Each cycle is 
replaced by a hyperbolic sink that lies close to the original cycle [3l|. For sufficiently small 
C^-perturbations, the existence of solutions that go several times around the cycles is ruled out. 

The fixed point hyperplane defined by Fix(( 72 )) = {{xj, X 2 , X 3 , X 4 ) G : X 3 = 0} divides 
in two flow-invariant connected components, preventing arbitrarily visits to both cycles in 
Trajectories whose initial condition lies outside the invariant subspaces will approach one of the 
cycles in positive time. Successive visits to both cycles require breaking this symmetry |34j . 

2.5. Breaking the Z 2 (( 72 ))-symmetry. From now on, we consider /o embedded in a generic 
one-parameter family of vector fields, breaking the ( 72 )-equivariance as follows: 

(P 6 ) The vector fields /a : —)• TS^ are a family of ( 7 i)-equivariant vector fields. 

Since the equilibria v and w lie on Fix(( 7 i)) and are hyperbolic, they persist for small A > 0 
and still satisfy Properties (F|2|) and (FjSj). Their invariant two-dimensional manifolds generically 
meet transversely along two trajectories. The generic bifurcations from a manifold are discussed 
by Chillingworth [T2|, under these conditions we assume: 

(P7) For A 7 ^ 0, the two-dimensional manifolds lT“(w) and lT®(v) intersect transversely at 
two trajectories that we will call primary connections. Together with the connections in 
(F|3|) this forms a Bykov heteroclinic network that we denote by 

The network consists of four copies of the simplest heteroclinic cycle between two saddle- 
foci of different Morse indices, where one heteroclinic connection is structurally stable and the 
other is not, a Bykov cycle. Property (F|7|) is natural since the heteroclinic connections of (F|Z|), 
as well as those of assertion ([3|) of Theorem [2] below, occur at least in symmetric pairs. In what 
follows, we describe the dynamics near a subnetwork consisting of the two primary connections 
together with one of the unstable connections. The results obtained concern each one of the two 
subnetworks of this type. 

For small A 7 ^ 0, the neighbourhood is still positively invariant and contains the network 
Since the closure of is compact and positively invariant it contains the w-limit sets of 
all its trajectories. The union of these limit sets is a maximal invariant set in V^. For /o this 
is the cycle S®, by Lemma (H whereas for symmetry-breaking perturbations of /o it contains S'*" 
but does not coincide with it. Our aim is to describe this set and its sudden appearance. 

We proceed to review the dynamics in a small tubular neighbourhood of the cycle. In order 
to do this we introduce some concepts. 

Let P C S'*" be one Bykov cycle involving v and w, and the connections [v —>■ w] and [w —>■ v] 
given by (F|3j) and (F|7|), respectively. Let V,W be disjoint neighbourhoods of the equilibria 
as above. Consider two local cross-sections of P at two points p and q in the connections [v —?■ w] 
and [w — > v], respectively, with p,q U W. Saturating the cross-sections by the flow, one 
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obtains two flow-invariant tubes joining V and W that contain the connections in their interior. 
We call the union of these tubes with V and W a tubular neighbourhood T of the Bykov cycle. 
More details will be provided in Section SI 

Definition 1. Let V,W C. be two disjoint neighbourhoods of v and w, respectively. A one¬ 
dimensional connection [w —v] that, after leaving W, enters and leaves both V and W precisely 
re G N times is called an re-pulse heteroclinic connection with respect to V and W. When there 
is no ambiguity, we omit the expression “with respect to V and W”. //re > 1 we call it a 
multi-pulse heteroclinic connection. If W“(w) and ly^(v) meet tangentially along [w —v], 
we say that the connection [w —)• v] is an re-pulse heteroclinic tangency, otherwise we call it a 
transverse re-pulse heteroclinic connection. 

The primary connections [w —> v] in of (F[7|) are transverse 0-pulse heteroclinic connec¬ 
tions. With these conventions we have: 

Theorem 2 ([M])- If ^ vector field fo satisfies (If^-(If^, then the following properties are sat¬ 
isfied generically by vector fields in an open neighbourhood of fo in the space of (ji)-equivariant 
vector fields of class on , r > 3: 

(1) there is a heteroclinic network S* consisting of four Bykov cycles involving two equilibria 
V and w, two 0-pulse heteroclinic connections [v —>■ w] and two 0-pulse heteroclinic 
connections [w ^ v]; 

(2) the only heteroclinic connections from v to w are those in the Bykov cycles and there 
are no homoclinic connections; 

(3) any tubular neighbourhood of a Bykov cycle T in S* contains infinitely many n-pulse 
heteroclinic connections [w v] for each re G N, that accumulate on the cycle; 

(4) for any tubular neighbourhood T, given a cross-section Sq C T at a point g ire [w —>■ v], 
there exist sets of points such that the dynamics of the first return to Sq is uniformly 
hyperbolic and conjugate to a full shift over a finite number of symbols. These sets 
accumulate on S* and the number of symbols coding the return map tends to infinity as 
we approach the network. 

Notice that assertion m of Theorem [2] implies the existence of a bigger network: beyond 
the original transverse connections [w ^ v], there exist infinitely many subsidiary heteroclinic 
connections turning around the original Bykov cycle. It also follows from m and 0 that 
any tubular neighbourhood 7~ of a Bykov cycle T in S* contains points not lying on T whose 
trajectories remain in T for all time. In contrast to the findings of [sn Ea [53], the suspended 
horseshoes of (jH) arise due to the presence of two saddle-foci together with transversality of 
invariant manifolds, and does not depend on any additional condition on the eigenvalues at the 
equilibria. 

A hyperbolic invariant set of a C'^-diffeomorphism has zero Lebesgue measure |H|- Never¬ 
theless, since the authors of |31| worked in the category, this set of horseshoes might have 
positive Lebesgue measure. Rodrigues m proved that this is not the case: 

Theorem 3 f[17]L Let T be a tubular neighbourhood of one of the Bykov cycles T of Theorem{^ 
Then in any cross-section Sq C T the set of initial conditions in Sq that do not leave T for all 
time, has zero Lebesgue measure. 

The shift dynamics does not trap most solutions in the neighbourhood of the cycle. In 
particular, none of the cycles is Lyapunov stable. 
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3 . Statement of results 

Heteroclinic cycles connecting saddle-foci with a transverse intersection of two-dimensional 
invariant manifolds imply the existence of hyperbolic suspended horseshoes. In our setting, 
when A varies close to zero, we expect the creation and the annihilation of these horseshoes. 
When the symmetry Z 2 (( 72 )) is broken, heteroclinic tangencies are reported in the next result. 

Theorem 4. In the set of families f\ of vector fields satisfying m-m there is a subset C, 
open in the topology, for which there is a sequence A* > 0 o/ real numbers, with limj_).oo A* = 0 
such that for X > A*, there are two 1-pulse heteroclinic connections for the flow ofx = f\{x), that 
collapse into a 1-pulse heteroclinic tangency at X = Xi and then disappear for A < A*. Moreover, 
the 1-pulse heteroclinic tangency approaches the original [w —)• v] connection when A* tends to 
zero. 

The explicit description of the open set C is given in Section [5l after establishing some notation 
for the proof. Although a tangency may be removed by a small smooth perturbation, the 
presence of tangencies is densely persistent, a phenomenon similar to the Cocoon bifurcations 
observed in the Michelson system [laiisj observed by Lau. 

Theorem 5. For a family f\ in the open set C of Theorem and for each parameter value 
Xi corresponding to a 1-pulse heteroclinic tangency, there is a sequence of parameter values Xij 
accumulating at Xi for which there is a 2-pulse heteroclinic tangency. This property is recursive 
in the sense that each n-pulse heteroclinic tangency is accumulated by {n + l)-pulse heteroclinic 
tangencies for nearby parameter values. 

Associated to the heteroclinic tangencies of Theorem [SJ suspended horseshoes disappear as A 
goes to zero. 

Theorem 6. For a family f\ in the open set C of Theorem there is a sequence of closed 
intervals = [cn,dn], with 0 < dn+i,Cn < dn and lim dn = 0, such that as X decreases in An, 

71^00 

a suspended horseshoe is destroyed. 

A similar result has been formulated by Newhouse m and Yorke and Alligood [53] for the 
case of two dimensional diffeomorphisms in the context of homoclinic bifurcations in dissipative 
dynamics, with no references to the equivariance. A more precise formulation of the result is 
given in Section [ 6 l Applying the results of [431 [53] to this family, we obtain: 

Corollary 7. For a family f\ in the open set C of Theorem^ and X in one of the intervals An 
of TheoremlE the flow of x = f\{x) undergoes infinitely many saddle-node and period doubling 
bifurcations. 

With the additional hypothesis that the first return map to a suitable cross-section contracts 
area, we get: 

Corollary 8. For a family f\ in the open set C of Theorem^ if the first return to a transverse 
section is area-contracting, then for parameters X in an open subset of An with sufficiently large 
n, infinitely many attracting periodic solutions coexist. 

In Section [ 6 | we also describe a setting where the additional hypothesis holds. When A de¬ 
creases, the Cantor set of points of the horseshoes that remain near the cycle is losing topological 
entropy, as the set loses hyperbolicity, a phenomenon similar to that described in [22|. For A ~ 0, 
return maps to appropriate domains close to the tangency are conjugate to Henon-like maps 
[I311291139]. As A —>■ 0, in V^, infinitely many wild attractors may coexist with suspended 
horseshoes that are being destroyed. 
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Figure 2. Cylindrical neighbourhoods of the saddle-foci w (a) and v (b). 

The complete description of the dynamics near these bifurcations is an unsolvable problem: 
arbitrarily small perturbations of any differential equation with a quadratic heteroclinic tangency 
may lead to the creation of new tangencies of higher order, and to the birth of quasi-stochastic 

attractors [201 [251 [231 [M!- 

4. Local geometry and transition maps 

We analyse the dynamics near the network by deriving local maps that approximate the 
dynamics near and between the two nodes in the network. In this section we establish the 
notation that will be used in the rest of the article and the expressions for the local maps. We 
start with appropriate coordinates near the two saddle-foci. 

4.1. Local coordinates. In order to describe the dynamics around the Bykov cycles of we 
use the local coordinates near the equilibria v and w introduced by Ovsyannikov and Shilnikov 
[42] . Without loss of generality we assume that a-v = Ow = 1 ■ 

In these coordinates, we consider cylindrical neighbourhoods V and W in of v and w, 
respectively, of radius p = e > 0 and height z = 2e — see Figure [2l After a linear rescaling of 
the variables, we may also assume that e = 1. Their boundaries consist of three components: 
the cylinder wall parametrised by x G R (mod 27r) and |y| < 1 with the usual cover 

{x,y) hA il,x,y) = {p,e,z) 

and two discs, the top and bottom of the cylinder. We take polar coverings of these disks 

(r,(p) hA (r,(p,±l) = {p,e,z) 

where 0 < r < 1 and p G R (mod 27r). The local stable manifold of v, IT^(v), corresponds to 
the circle parametrised hy y = 0. In V we use the following terminology suggested in Figure [21 

• In(v), the cylinder wall of V, consisting of points that go inside V in positive time; 

• Out(v), the top and bottom of V, consisting of points that go outside V in positive time. 
We denote by In''"(v) the upper part of the cylinder, parametrised by (x,y), y G [0,1] and by 
In“(v) its lower part. 

The cross-sections obtained for the linearisation around w are dual to these. The set ly^(w) 
is the z-axis intersecting the top and bottom of the cylinder W at the origin of its coordinates. 
The set IT“(w) is parametrised hy z = 0, and we use: 
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• In(w), the top and bottom of W, consisting of points that go inside W in positive time; 

• Out(w), the cylinder wall of W, consisting of points that go inside W in negative time, 
with Out''"(w) denoting its upper part, parametrised by (x,y), y G [0,1] and Out“(w) 
its lower part. 

We will denote by portion of W‘^{'w) that goes from w up to In(v) not intersecting 

the interior of V and by portion of W^(v) outside W that goes directly from Out(w) 

into V. The flow is transverse to these cross-sections and the boundaries of V and of W may be 
written as the closure of In(v) U Out(v) and In(w) U Out(w), respectively. 


4.2. Local maps near the saddle-foci. Following |14l 142]. the trajectory of a point (x,y) 

with y > 0 in In(v), leaves V at Out(v) at 

(4.2) 


4>v(x,y) = { +S'i(x,y;A),-^ + x + S’ 2 (x, y; A) ) = (r, 0) 


where 


A 


where and S 2 are smooth functions which depend on A and satisfy; 

gk+l+r> 


(4.3) 


dx^dx^dX^ 


Si{x,y;X) 


<cy 


.S^+cr—l 


and C and a are positive constants and k, I, m are non-negative integers. Similarly, a point (r, 4>) 

in In(w)\lTj®^(w), leaves W at Out(w) at 

(4.4) 


^w(i’,¥5) = ( -h¥?-Fi?i(r,v2;A),/"' -h i?2(r, A) ) = (x,y) 


where 


<5w = ^ > 1 

Xj-oi 


where Ri and R 2 satisfy a condition similar to (14.21) . The terms Si, S 2 , R 2 correspond to 
asymptotically small terms that vanish when y and r go to zero. A better estimate under a 
stronger eigenvalue condition has been obtained in [271 Prop. 2.4]. 


4.3. Transition map along the connection [v —wj. Points in Out(v) near 
mapped into In(w) in a flow-box along the each one of the connections [v —>■ w]. Without loss 
of generality, we will assume that the transition Tv-j.w : Out(v) —>■ ln(w) does not depend on A 
and is modelled by the identity, which is compatible with hypothesis (Fl5]). Using a more general 
form for Tv-s-w would complicate the calculations without any change in the final results. 

The coordinates on V and W are chosen to have [v —>■ w] connecting points with 2 > 0 in U 
to points with z > 0 in W. We will denote by t] the map r] = o 'I'v^w ° ‘l*v Up to higher 
order terms, from ()4.2p and ()4.4p . its expression in local coordinates, for y > 0, is 

(4.5) r/(x,y) = ("x — iP In y,y'^) with 5 = (5v5w > 1 and K = ^ > 0. 

The choice of as the identity reflects the fact that the node have the same chirality. 

In the case where the nodes have different chirality, them map Tv-i,w reverses orientation. This 
affects the form of the map y and any subsequent results that use it. 


4.4. Geometry of the transition map. Consider a cylinder C parametrised by a covering 
(0,y) G R X [—1,1], where 6 is periodic. A helix on the cylinder C accumulating on the circle 
y = 0 is a curve on C without self-intersections, that is the image, by the parametrisation {9, y), 
of a continuous map H : {a,b) — )• R x [—1,1], H{s) = {H 0 {s), Hy{s)), satisfying: 

lim Hq{s) = lim Hq{s) = ±00 lim Hy{s) = lim Hy[s) = 0 

s—>-a+ s^b~ s—>-a+ s^b~ 
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and such that there are a < 6 G icL,b) for which both Hg{s) and Hy{s) are monotonic in each 
of the intervals (a, a) and (b,b). It follows from the assumptions on the function Hg that it has 
either a global minimum or a global maximum, since lims_j,„+ Hg{s) = lim^^^- Hg(s). At the 
corresponding point, the projection of the helix into the circle y = 0 is singular, a fold point 
of the helix. Similarly, the function Hy always has a global maximum, that will be called the 
maximum height of the helix. See Figure [3l 




Figure 3. A helix is defined on a covering of the cylinder by a smooth curve 
{Hg{s), Hy{s)) that turns around the cylinder infinitely many times as its height Hy 
tends do zero. It always contains a fold point, indicated here by a diamond, and a point 
of maximum height, shown as a black dot. 


Lemma 9. Consider a curve on one of the cylinder walls In{v) or Out('w), parametrised by the 
graph of a smooth function h : [o, 6] —R, where b — a < 2 'k with h{a) = h{b) = 0, h'{a) > 0, 
h'{b) < 0 and h{x) > 0 for all x € {a,b). Let M he the global maximum value of h, attained at 
a point xm S Then, for the transition maps defined above, we have: 

( 1 ) if the curve lies in In{y) then it is mapped by rj = o'I'v^.w o‘hv into a helix on Out{w) 

accumulating on the circle Out{w) n its maximum height is and it has a 

fold point at y(x*,/i(x*)) for some x* G {a,XM); 

( 2 ) if the curve lies in Out{w) then it is mapped hy rj~^ into a helix on In{'v) accumulating 

on the circle In{v) n its maximum height is and it has a fold point at 

ri~^{x^, /i(x*)) for some x* G (xm, b). 

This result depends strongly on the form of rj, and therefore, on the chirality of the nodes. 
If the nodes had different chirality, the graph of h would no longer be mapped into a helix; the 
curve rj (x, h{x)) would instead have a vertical tangent at infinitely many points, see [M] . 

Proof. The graph of h defines a curve on In(v) without self-intersections. Since rj is the transition 
map of a differential equation, hence a diffeomorphism, this curve is mapped by r] into a curve 
H{x) = rj{x,h{x)) = {Hg{x),Hy{x)) in Out(w) without self-intersections. Using the expression 
(j4.5p for y, we get 

(4.6) R(x) = ^x — iFln/i(x), (h(x))^^ and H'{x) = 

From this expression it is immediate that 

lim Hg{x) = lim Hg{x) = -boo and lim Hy{x) = lim Hy{x) = 0. 

x—>-a+ x—>^b~ ai—>-a+ x—>b~ 


f _ Kh\x) Jh{x)f h'{x) \ 
h{x) ’ h{x) j ■ 


Also, 


Hy{XM) = > {h{x)f = Hy{x) 
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for all X G {a,b), so the curve lies below the level y = in Out(w). Since h'{b) < 0, there 
is an interval {b,b) where h'(x) < 0 and hence H'g{x) > 0 and Hy{x) < 0. Similarly, h'{x) > 0 
on some interval {a, a), where Hy{x) > 0. For the sign of Hq{x), note that Kh'{a) > 0 = h[a) 
and Kh'{xM) = 0 < h{xM)- Thus, there is a point x* G {a,XM) where /i(x*) = Kh'{x^,) and 
x) changes sign, this is a fold point. If x* is the minimum value of x for which this happens, 
then locally the helix lies to the right of this point. This proves assertion ([T]). The proof of 
assertion Q is similar, using the expression 


(4.7) 


K 


r] {x,y) = [x+ — liiy,y 


1/5 


In this case we get 


lim Hg{x) = lim Hq{x) = —oo 

x—>-a+ x^b~ 

and if x* is the largest value of x for which the helix has a fold point, then locally the helix lies 
to the left of r/“^(x*,/i(x*)). □ 


Note that if, instead of the graph of a smooth function, we consider continuous, positive and 
piecewise smooth curve a{s) without self-intersections from a(0) = (o, 0) to a(l) = (6,0), we 
can apply similar arguments to show that both r/(a) and rj~^{a) are helices. 


4.5. Geometry of the invariant manifolds. There is also a well defined transition map 

■ Out(w) —^ In(v) 

that depends on the Z2(72)-symmetry breaking parameter A, where is the identity map. 

We will denote by R\ the map o rj, where it is well defined. When there is no risk of 

ambiguity, we omit the superscript A. In this section we investigate the effect of on the 

two-dimensional invariant manifolds of v and w for A 7^ 0, under the assumption (F[7]). 

For this, let f\ be an unfolding of /o satisfying (F[I])-(F1Z1). For A 7^ 0, we introduce the 
notation, see Figured) 

• (P.^, 0) and (P^, 0) with 0 < P^ < P^ < 27r are the coordinates of the two points where 
the connections [w —>■ v] of Property (F)?]) meet Out(w); 

• (Pv,0) and (Pv,0) with 0 < Pj: < P^ < 2?! are the coordinates of the two points where 
[w v] meets In(v); 

• (Pw,0) and (Pv,0) are on the same trajectory for each j = 1,2. 

By (FIZD, for A 7^ 0, the manifolds iy“(w) and IF*(v) intersect transversely along the primary 
connections. For A close to zero, we are assuming that ITioc,(v) intersects the wall Out(w) of 
the cylinder VF on a closed curve as in Figured) It corresponds to the expected unfolding from 
the coincidence of the manifolds ly^(v) and iy“(w) at /q. Similarly, ITioc(w) intersects the wall 
In(v) of the cylinder V on a closed curve. For small A > 0, these curves can be seen as graphs 
of smooth 27r-periodic functions, for which we make the following conventions: 

• IFi*^(v) n Out(w) is the graph of y = /iv(x. A), with /iv(Pw, A) = 0, j = 1, 2; 

• IF[“^(w) n In(v) is the graph of y = /iw(a:. A), with /iw(Pv, A) = 0, j = 1,2; 

• /iv(x,0) = 0 and hv,{x,0) = 0 

• for A > 0, we have /i(,(P^, A) > 0, hence A) < 0 and A) < 0, h'^{P^, A) > 

0. 

The two points (P^,0) and (P^,0) divide the closed curve y = /iv(x,A) in two components, 
corresponding to different signs of the second coordinate. With the conventions above, we get 
/iv(x,A) > 0 for X G (P^,P^). Then the region W~ in Out(w) delimited by ITioc(^) 
IFj“c(w) between P^ and P^ gets mapped by into In~(v), while all other points in 
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Figure 4. For A close to zero, VF®(v) intersects the wall Out(w) of the cylinder W 
on a closed curve, given in local coordinates as the graph of a periodic function. Simi¬ 
larly, fF“(w) meets In(v) on a closed curve — this is the expected unfolding from the 
coincidence of the invariant manifolds at A = 0. 

Out“'"(w) are mapped into In^(v). We denote by W~^ the latter set, of points in Out(w) with 
0 < y < 1 and y > hv(x,\) for x G The maximum value of hv{x,X) is attained at 

some point 

(x,y) = (xv(A),Mv(A)) with P^ < x^iX) < P^. 

Finally, let M^{X) be the maximum of h^{x,X), attained at a point a:w(A) G In 

order to simplify the writting, the following analysis is focused on the two Bykov cycles whose 
connection [v — >■ w] lies in the subspace defined by X 3 > 0. The same results, with minimal 
adaptations, hold for cycles obtained from the other connection. With this notation, we have: 

Proposition 10. Let f\ be family of vector fields satisfying m-m- For A / 0 sufficiently 
small, the portion of n /n(v) that lies in is mapped by rj into a helix in Out{w) 

accumulating on IF;“^(w). If Mw(A) is the maximum height of n /n“''(v), then the 

maximum height of the helix is Mw(A)‘^. For each A > 0 there is a fold point in the helix that, 
as X tends to zero, turns around the cylinder wall Out{'w) infinitely many times. 

Proof. That rj maps IT[“^(w) nln^(v) into a helix, and the statement about its maximum height 
follow directly by applying assertion ([1]) of Lemma [9] to /iw For the fold point in the helix, let 
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x*(A) be its first coordinate. From the expression (14.51) of rj it follows that 

x*(A) = x\ — K\n.h^{x\) 

for some x\ G (Py,Xw(A)) and with h^{x\, A) < hy,{x^{X), A) = M^{X) and hence, 

-K\nh^{xx,X) > -KlnMw(A). 

Since fx unfolds /o, then hm;^_j,o iWw(A) = 0, hence lim;,^o —Kin h^{xx, A) = oo and therefore, 
the fold point turns around the cylinder Out(w) infinitely many times. □ 

The Hypothesis (Fl5]) about chirality of the nodes is essential for Lemma[9]and PropositionllOl 
If we had taken the rotation in W with the opposite orientation to that in V, the rotations would 
cancel out and the curve H{x) dehned in (14.6p would no longer be a helix, since the angular 
coordinate Hg would not be monotonic. 


5. Heteroclinic Tangencies 

Using the notation and results of Section 0] we can now discuss the tangencies of the invariant 
manifolds and prove Theorems 0] and El As remarked in Section 14.51 since fx unfolds /o, then 
the maximum heights, Mv(A) of Out(w), and Mw(A) of Wi“^(w) n In(v), satisfy: 

lim Mv(A) = lim M^fA) = 0. 

A^-o \->-o 

We make the additional assumption that (Mw(A))‘^ tends to zero faster than Mv(A). This 
condition dehnes the open set C of unfoldings fx that we need for the statement of Theorem 01 
The subset C of Theorem 01 is the set of families fx of vector fields satisfying (F[I])-(Fl7]) for 
which there is a value A* > 0 such that for 0 < A < A* we have (Mw(A))^ < Mv(A). Then C is 
an open subset, in the topology, of the set of families fx of vector helds satisfying (F[I])-(Fl71). 


Out(w) 


ln(v) 



X 


Figure 5. Left: when A decreases, the fold point of the helix a^{x) G Out(w) moves 
to the right and for A = A^ it is tangent to Wi®j.(v) creating a 1-pulse tangency. Right: 
^'[w^v] niaps the helix close to Wi“^(w) n In(v), creating several curves that 

satisfy the hypotheses of Lemma 01 (H} ■ These curves are again mapped by rj into helices 
in Out(w) creating 2-pulse heteroclinic tangencies. 
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5.1. Proof of Theorem [4l Suppose f\ G C. By Proposition 1101 the curve 

a^{x) = rj {x,hv,{x, X), X) = (x), 02 (x)^ , x G (mod 27r) 

is a helix in Out'^('w) and has at least one fold point at x = x*(A). The second coordinate of 
the helix satisfies 0 < < XI^{XY for all x G (P^,P^) (mod 27r) and all positive A < A*. 

Since f\ G C, then a^ix) < Mv(A) for all x and all positive A < A*. 

Moreover, since the fold point a^(x*(A)) turns around Out(w) infinitely many times as A goes 
to zero, given any Aq < A* there exists a positive value Xn < Aq such that a^(x*(Ai?)) lies in 
W~^, the region in Out(-w) between lT[®^(v) and lT[“(,(w) that gets mapped into the upper part 
of In(v). Since the second coordinate of the fold point is less than the maximum of /iv, there 
is a positive value Xl < Xr such that a^{x^,{XL)) lies in W~, the region in Out(w) that gets 
mapped into the lower part of In(v), whose boundary contains the graph of h^. Therefore, the 
curve a'^(x*(A)) is tangent to the graph of /iv(x, A) at some point a(x*(Ai)) with Ai G (Xl, Xr). 

We have thus shown that given Aq > 0, there is some positive Ai < Aq, for which the image of 
the curve lT[“^(w) n In(v) by rj is tangent to Wiq^(v) n Out(w), creating a 1-pulse heteroclinic 
tangency. Two transverse 1-pulse heteroclinic connections exist for A > Ai close to Ai. These 
connections come together at the tangency and disappear. 

By Proposition [ini as A goes to zero, the fold point a^(x*(A)) turns around the cylinder 
Out(w) infinitely many times, thus going in and out of W~. Each times it crosses the boundary, a 
new tangency occurs. Repeating the argument above yields the sequence A* of parameter values 
for which there is a 1-pulse heteroclinic tangency and this completes the proof of the main 
statement of Theorem 01 

On the other hand, as A goes to zero, the maximum height of the helix, (Mw(A))'^ also tends to 
zero. This implies that the second coordinate of the points a'’'(x*(Aj)) G Out(w) where there is 
a 1-pulse heteroclinic tangency tends to zero as i goes to infinity. This shows that the tangency 
approaches the two-dimensional connection [w —>■ v] that exists for A = 0. □ 

A crucial fact in the proof of Theorem 0] is that W^^('w) n In(v) is the graph of a function 
with a single maximum, and that this maximum goes to 0 as A goes to 0. This follows because 
the family fx unfolds the more symmetric vector field /q. We could make the assumption on the 
invariant manifold directly, but it is the context of symmetry-breaking that makes it natural. 

The construction in the proof of Theorem 0] may be extended to obtain multipulse tangencies, 
as follows: 

5.2. Proof of Theorem [51 Look at W[“^(w) nln(v), the graph of /iw(a^, A). Since h'^{P^,X) < 
0, there is an interval [x, P^] C [xw(A), where the map is monotonically decreasing. Note 
that the hypothesis of different chirality is essential here. Therefore, we may define infinitely 
many intervals where a'^(x) = rj (x, h^(x, A)) lies in W~^. More precisely, we have two sequences, 
(aj) and (hj) in [x,Pj] such that: 

• ttj < bj < Oj+i with limj^oo Oj = Pj; 

• a^(aj) and a^(hj) G n Ouf'"(w); 

• if X G (aj,bj) then a'^(x) G W^; 

• the curves {[aj, bj]) accumulate uniformly on ^[“(.(w) n Out(w) as j 00. 

Hence, each one of the curves ([a^, bj]) is mapped by 'kw-^v into the graph of a function h in 
In(v) satisfying the conditions of Lemma [9l and hence each one of these curves is mapped by rj 
into a helix ^j(x), x G (uj, bj). 

Let Aj be a parameter value for which x = fx^ (x) has a 1-pulse heteroclinic tangency as stated 
in Theorem 01 As j —>■ -|-oo, the helices ^j(x) accumulate on the helix of Theorem 01 as drawn in 
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Out{w) 


1 


^j(x) 




a(x) 


pi p2 

\ 

w' (v) n Out(w) 

loc 

W'" (w)nOut(w) 
loc 


Figure 6 . Curves in the proof of Theorem[5] for X = Xi the curve a^{x) = i?(Wj“(,(w)n 
In'''(v)) (solid black curve) is tangent to (gray curve) in Out(w). The curves 

^j{x) (dotted) accumulate on a^(x). Very small changes in A make them tangent to 
IFjqc(v) creating 2-pulse heteroclinic tangencies. 


Figure[6l hence the fold point of ^j(x) is arbitrarily close to the fold point of r/(x, hv,{x, A)). The 
arguments in the proof of Theorem [3] show that a small change in the parameter A makes the new 
helix tangent to 1 Fjqj,(v) as in Figure [5j For each j this creates a 2-pulse heteroclinic tangency 
at A = Xij. Since the the helices Cji^) accumulate on a^*(x), it follows that lim^gN ^ij = A*. 

Finally, the argument may be applied recursively to show that each n-pulse heteroclinic 
tangency is accumulated by (n -|- l)-pulse heteroclinic tangencies for nearby parameter values. 

□ 

If A* G R is such that the flow of x = fx* {x) has a heteroclinic tangency, when A varies near 
A*, we find the creation and the destruction of horseshoes that is accompanied by Newhouse 
phenomena [221 [53]. In terms of numerics, we know very little about the geometry of these at¬ 
tractors, we also do not know the size and the shape of their basins of attraction. Basin boundary 
metamorphoses and explosions of chaotic saddles as those described in [35| are expected. 

6 . Bifurcating dynamics 

We discuss here the geometric constructions that determine the global dynamics near a Bykov 
cycle, in order to prove Theorem [6j For this we need some preliminary definitions and more 
information on the geometry of the transition maps. First, we adapt the definition of horizontal 
strip in |26] to serve our purposes: for r > 0 sufficiently small, in the local coordinates of the 
walls of the cylinders V and W, consider the rectangles: 

= [-Pv -T,^v X [0) 1] C In(v) and <Sw = [P^ - r, -b r] x [0,1] C Out(w) 

with the conventions —tt <Pv“''‘<'Pv+'^^^ —vr < — t < + r < tt. 

A horizontal strip in <Sv is a subset of In'’'(v) of the form 

"P = {(x,y) €/n+(v) : x G - r, -b r] and y G [tti(x), U2 (x)]} , 

where ui,tt2 : [Pv — t, Pj -b r] (0,1] are smooth maps such that ui{x) < U 2 {x) for all x G 
— r, Py -b r]. The graphs of the Uj are called the horizontal boundaries of P and the segments 
(p2 -T,y), ui(p2 -t) <y < U2(Pv - r) and (P^ + T,y), with ui(Pj + t) < y < U2{P^ + r) 
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are its vertical boundaries. Horizontal and vertical boundaries intersect at four vertices. The 
maximum height and the minimum height of T-L are, respectively 

max U2(x) min ui(x). 

xG[P2_T,P^+r] xG[P3-T,Pl+r] 

Analogously we may define a horizontal strip in 5w C Ouf^( w). 

A horseshoe strip in <Sv is a subset of In“’“(v) of the form 

{(x,y) e In+(v) : xG [02,62], y G [oi(x),02(x)]} , 

where 

• [oi, 61] C [02, 62] C [Pv -T,Pv +'ry, 

• oi(ai) = 01(61) = 0; 

• 1^2(02) = 1^2(62) = 0. 

The boundary of T-L consists of the graph of U2{x), x G [02,62], the graph of ui{x), x G [oi,6i], 
together with the two pieces of ITjq^(v) D In(v) with x G [02,01] and x G [61,62]. 

6 . 1 . Horseshoe strips in 5 v. We are interested in the dynamics of points whose trajectories 
start in and return to In^(v) arriving at S^. 


In{v) 



Figure 7. For A > 0 sufficiently small, the set 77 “^(5w) H has infinitely many 
connected components (gray region), all of which, except maybe for the top ones, define 
horizontal strips in <Sv accumulating on W;®^(v) n In(v) (gray curve). 

Lemma 11. The set ?7“^(5w) H 5v has infinitely many connected components all of which are 
horizontal strips in accumulating on n/n(v), except maybe for a finite number. The 

horizontal boundaries of these strips are graphs of monotonically increasing functions of x. 

Proof. The boundary of 5w consists of the following: 

(1) a piece of lT[“^(w) n Out(w) parametrised by y = 0, x G [P^ — r, P^ + r], where r]~^ is 
not defined; 

(2) the horizontal segment (x, 1) with x G [P^ — T,Pl+T-]-, 

(3) two vertical segments (P^ — r, y) and (P^ + r, y) with y G (0,1). 

Together, the components m and dSD form a continuous curve that, by the arguments of 
Lemma [9] (l2|), is mapped by y into a helix on In”'"(v), accumulating on lT[®^(v) n In(v). As 
the helix approaches lT[®(,(v), it crosses the vertical boundaries of <Sv infinitely many times. 
The interior of 5w is mapped into the space between consecutive crossings, intersecting <Sv in 
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horizontal strips, as shown in Figure [71 From the expression (|4.7p of it also follows that 
the vertical boundaries of <Sw are mapped into graphs of monotonically increasing functions of 

X. □ 

Here, as in the next two Lemmas, we are using the form of the map rj through Lemma [9l 
hence the result depends strongly on the chirality hypothesis (P|5]). 

Denote by Hn the strip that attains its maximum height at the vertex + r, hn) with 


( 6 . 8 ) 


U _ {Pi-P^+2T-2nn)/K 
icn — c , 


then lim„_,.oo hn = 0, hence the strips Hn accumulate on VFiqc(v) H In(v). The minimum height 
of Hn is given by 




(6.9) 

and is attained at the vertex (P^ — r, m„). Moreover n < m implies that Hn lies above Hr^ 


ln(v) 



a 

b 


Figure 8 . Horseshoe strips: for A > 0 sufficiently small, the set R\{r] ^(5w) n5v) C 
In"^ (v) defines a countably infinite number of horseshoe strips in S-v , that accumulate on 
VFi“,(w)nln(v). 


Lemma 12. Let Hn be one of the horizontal strips in t/“^(5w) n<Sv. Then r](Hn) is a horizontal 
strip in 5w. The strips rj{Hn) accumulate on W^^{w) n /n(v) as n ^ oo and the maximum 
height of r]{Hn) is hf^, where hn is the maximum height ofHn- 

Proof. The boundary of consists of a piece of VLiq^(v) plus a curve formed by three segments, 
two of which are vertical and a horizontal one. From the arguments of Lemma [9] ([T]), it follows 
that the part of the boundary of not contained in Wjq^(v) is mapped by g into a helix. 

Consider now the effect of g on the boundary of Hn- Each horizontal boundary gets mapped 
into a piece of one of the vertical boundaries of S-w- The vertical boundaries of Hn are contained 
in those of S-v and hence are mapped into two pieces of a helix, that will form the horizontal 
boundaries of the strip g{Hn), that may be written as graphs of decreasing functions of x. 
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A shown after Lemma m the maximum height of Tin tends to zero, hence the strips 'r]{'Hn) 
have the same property. The maximum height of r]{'Hn) is hf^, attained at the point (P.^ — r, hi). 

□ 

Lemma 13. For each A > 0 sufficiently small, there exists n-o(A) such that for all n > no the 
image R\{'Hn) of the horizontal strips in 77“^(5w) n <Sv intersects in a horseshoe strip. The 
strips R\{'Hn) accumulate on n /n(v) and, when n 00 , their maximum height tends 

to Mw(A), the maximum height of n /n(v). 

Proof. The curve ILioc(v) H Out(w) is the graph of the function hv{x,X) that is positive for x 
outside the interval [P^, P^] (mod 27r). In particular, h^{P^ — r. A) > 0 and hv{P^ + l. A) >0 
for small r > 0 . Therefore there is a piece of the vertical boundary of 5w that lies below 
^ioc(^) Out(w), consisting of the two segments (P^ — T,y) with 0 < y < h^{P‘^ — r. A) and 
(P^ + T, y) with 0 < y < h^{Pi, + r. A). For small A > 0 , these segments are mapped by 
inside In~(v). 

Let no be such that the maximum height of rjiRno) than the minimum of 

hv{P^ — T, A) > 0 and h^{P‘^ + t, A) > 0. 

Then for any n > no the vertical sides of ri{PLn) are mapped by inside In“(v). The 

horizontal boundaries of rjifHn) go across 5w, so writing them them as graphs of ui{x) < U 2 {x), 
there is an interval where the second coordinate of Uj{x)) is more than Mw(A) > 0, the 

maximum height of IT[“(,('w) fl In'’'(v). Since the second coordinate of Uj(x)) changes 

sign twice, then it equals zero at two points, hence R\{'Hri) n 5v is a horseshoe strip. 

We have shown that the maximum height of giTln) tends to zero as n —)• 00, hence the 
maximum height of R\{Rn) = ^w^.v(^(^n)) tends to Mw(A). □ 

6.2. Regular Intersections of Strips. We now discuss the global dynamics near the Bykov 
cycle. The structure of the non-wandering set near the network depends on the geometric 
properties of the intersection of Tin and R\{'Hri). 

Let A be a horseshoe strip and P be a horizontal strip in 5v We say that A and B intersect 
regularly if A n P 7^ 0 and each one of the horizontal boundaries of A goes across each one of 
the horizontal boundaries of P. Intersections that are neither empty nor regular, will be called 
irregular. 

If the horseshoe strip A and the horizontal strip P intersect regularly, then A n P has at least 
two connected components, see Figure [8l In this and the next subsection, we will find that the 
horizontal strips Fin across Pv may intersect R\{Rm) in the three ways: empty, regular and 
irregular, but there is an ordering for the type of intersection, as shown in Figure [9l 

Lemma 14. For any given fixed A > 0 sufficiently small, there exists N(\) € N such that for 
all m,n > N(X), the horseshoe strips Rx{FLm) in Pv intersect each one of the horizontal strips 
Rn regularly. 

Proof. From Lemma [T3] we obtain no(A) such that all R\{FLm) with m > no(A) are horseshoe 
strips and their lower horizontal boundary has maximum height bigger than the maximum height 
A/w(A) of IFi“^(w) n In(v). 

On the other hand, since the strips Fin accumulate uniformly on IFj®^(v) n In(v), with their 
maximum height hn tending to zero, then there exists ni(A) such that hn < iWw(A) for all n > 
ni(A). Note that the hn do not depend on A. Therefore, for m,n > N{X) = max{no(A),ni(A)} 
both horizontal boundaries of R\{FLm) go across the two horizontal boundaries of Fin- D 

The constructions of this section also hold for the backwards return map R~^ with analogues 
to Lemmas and M 
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Figure 9. Strips in Proposition [TS] (a) A = A3; (b) A = A2; (c) A = Ai. The position 
of the horizontal strips 'Ha and Hb does not depend on A. The maximum height of the 
horseshoe strips R\{Ha) and R\{Hb) decreases when A decreases, and the suspended 
horseshoe in R\[Ha) H Ha is destroyed. 


Generically, for n > N{X) each horizontal strip Hn intersects R\{Hn) in two connected 
components. Thus the dynamics of points whose trajectories always return to In(v) in Hn may be 
coded by a full shift on two symbols, that describe which component is visited by the trajectory 
on each return to Hn- Similarly, trajectories that return to inside Hn U • • • U Hn+k niay be 
coded by a full shift on 2k symbols. As fc —>• 00, the strips Hn+k approach VFi®^(v) n In(v) and 
the number of symbols tends to infinity. We have recovered the horseshoe dynamics described 
in assertion ([3|) of Theorem [2j 

The regular intersection of Lemma [TT] implies the existence of an i?y-invariant subset in 
17 “^ (5w) n5v, the Cantor set of initial conditions; 

^=0 U {Rx'{'»-m)nHn), 

jGZ m,n>N{X) 

where the return map to r]~^ (iSw)n5v is well defined in forward and backward time, for arbitrarily 
large times. We have shown here that the map R\ restricted to this set is semi-conjugate to a full 
shift over a countable alphabet. Results of [3l[3l] show that the first return map is hyperbolic 
in each horizontal strip, implying the full conjugacy to a shift. The time of return of points in 
Hj tends to + 00 , as j —>■ + 00 . The set A depends strongly on the parameter A, in the next 
subsection we discuss its bifurcations when A decreases to zero. 

6.3. Irregular intersections of strips. The horizontal strips Hn that comprise ? 7 “^(iSw) n5v 
do not depend on the bifurcation parameter A, as shown in Lemmas [TT] and [T 2 j This is in 
contrast with the strong dependence on A shown by the first return of these points to at the 
horseshoe strips R\{Hn)- In particular, the values of reo(A) fLemma fT3]l and N{X) (LemmadH) 
vary with the choice of A. For a small fixed A > 0 and for m,n > N(X) we have shown that Hn 
and R\{Hm) intersect regularly. 

The next result describes the bifurcations of these sets when A decreases. These global 
bifurcations have been described by Pahs and Takens in [l3] in a different context, where the 
horseshoe strips are translated down as a parameter varies. In our case, when A goes to zero the 
horseshoe strips are flattened into the common invariant two-dimensional manifoldsof v and w. 

Proposition 15. Given A 3 > 0 sufficiently small, there exist Ai < A 2 < A 3 € a horizontal 
strip Ha across C In'^ (v) and bo > a such that for any b > the horizontal strips Ha and 
Hb satisfy: 

(1) for A = A 3 the sets Hi and Rx^{Hj) intersect regularly for i,j € {a, b}; 

(2) for A = A 2 the intersection Ha H R\ 2 {Ha) is irregular; 

(3) for A = Ai the sets Ha and R\j^{Ha) do not intersect at al; 

(4) for A = Al and A = A 2 the set R\{Hb) intersects both Hb and Ha regularly. 
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Proof. As we have remarked before, the horizontal strips Tin do not depend on the parameter 
A. In particular, they are well defined for A = 0. Their image by the return map Rq is no 
longer a horseshoe strip, it is a horizontal strip across Sv Since for A = 0 the map is the 

identity (see 14.51) and the network is asymptotically stable, it follows that the maximum height 
of Ro{Rn) is no more than that of r](R.n)- 

From Lemma [12] and the expressions (|6.8p and (|6.9p for the maximum, hn, and minimum m^, 
height of Rn, we get that the maximum height /i^ of r]{Rn) is less than if 

Ln = 5(Pv - Pi) + 2t(<5 - 1) - 2 n 7 r (,5 - 1) < - P^. 

Since lim^^oo Ln = —oo, there is n 2 such that Ro{Rn) H Rn = 0, for every n > n 2 . Therefore, 
since Rx is continuous on A, for each n > n 2 there is A*(n) > 0 such that the maximum height 
of R\(n){Lin) is less than that of Rn- 

Let a = max{ni, A^(A 3 )} for N^X^) from Lemma fTTl Then assertion ([T|) is true for this Ra 
and for any Rn with n > a, by Lemma [TT| We obtain assertion Q by taking Ai = A*(a). 
Assertion ([2|) follows from the continuous dependence of R\ on A. Assertion (Uj) holds for any 
Rb with b > N{Xi) = bo by Lemma [TH □ 

6.4. Proof of Theorem [6] In particular, Proposition [15] implies that there exists d £ [A 2 , A 3 ) 
such that at A = d, the lower horizontal boundary of R\{Ra) is tangent to the upper horizontal 
boundary of Ra- Analogously, there exists c £ (Ai, A 2 ] such that at A = c, the upper horizontal 
boundary of Rx(Ra) is tangent to the lower horizontal boundary of Ra- There are infinitely 
many values of A in [c, d] for which the map Rx has a homoclinic tangency associated to a 
periodic point see [451IM] . 

The rigorous formulation of Theorem [H] consists of Proposition [TJ] , with Ai = [c, d] as in the 
remarks above. Since Proposition [TJ] holds for any A 3 > 0, it may be applied again with A 3 
replaced by Ai, and the argument may be repeated recursively to obtain a sequence of disjoint 
intervals A„. □ 

6.5. Proof of Corollary [3. The A dependence of the position of the horseshoe strips is not a 
translation in our case, but after Proposition [TJ] the constructions of Yorke and Alligood [Sl| 
and of Palis and Takens [l3] can be carried over. Hence, when the parameter A varies between 
two consecutive regular intersections of strips, the bifurcations of the first return map can be 
described as the one-dimensional parabola-type map in [l3| . Following |43l IM] , the bifurcations 
for A £ [Ai, A 3 ] are: 

• for X > d: the restriction of the map Rx^ to the non-wandering set on Ra is conjugate 
to the Bernoulli shift of two symbols and it no longer bifurcates as A increases. 

• at A £ (c, d): a fixed point with multiplier equal to ±1 appears at a tangency of the 
horizontal boundaries. It undergoes a period-doubling bifurcation. A cascade of period¬ 
doubling bifurcations leads to chaotic dynamics which alternates with stability windows 
and the bifurcations stop at A = d; 

• at A < c: trajectories with initial conditions in Ra approach the network and might be 
attracted to another basic set. 

□ 

6 . 6 . Proof of Corollary [ 8 ]. If the first return map Rx is area-contracting, then the fixed point 
that appears for A £ (c, d) is attracting for the parameter in an open interval. 

This attracting fixed point bifurcates to a sink of period 2 at a bifurcation parameter A > c 
close to c. This stable orbit undergoes a second flip bifurcation, yielding an orbit of period 4. 
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This process continues to an accumulation point in parameter space at which attracting orbits 
of period 2^ exist, for all A: G N. This completes the proof of Corollary [8] . □ 

Finally, we describe a setting in which the fixed point that appears for A G (c, d) can be shown 
to be attracting. Recall that the set Ir(v) is the graph of y = A). Suppose the 

transition map : Out(w) —In(v) is given by = {x,y + A)), which is 

consistent with Section IT31 Since 


J j and Dr]{x,y) = 

then 

det DRxix, y) = det (y(x, y)) ■ det Dri{x, y) = 5y^~^. 

For sufficiently small y (in T-Ln with sufficiently large n) this is less than 1, and hence R\ is con¬ 
tracting. However, if the first coordinate of y) depends on y, then det v)) 

will contain terms that depend on x -|- R In y and a more careful analysis will be required. 
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